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SYNTHESIS  OF  MCJ  CODES 


INTRODUCTION 

It  has  become  increasingly  important  to  employ  an  effective  mechanism  for  error  con¬ 
trol  in  the  initial  design  or  redesign  of  communication  systems,  both  for  channel  protection 
and  for  end-to-end  use.  The  most  efficient  codes  for  these  applications  are  the  maximum 
distance  separable  (MDS)  codes,  and  the  intent  here  is  to  present  constructive  techniques 
for  the  implementation  of  Massey -Costello-Justesen  (MCJ)  codes,  an  important  class  of  MDS 
codes. 


SYNTHESIS 

The  purpose  of  this  study  was  to  improve  our  understanding  of  constructive  techniques 
for  the  engineering  implementation  of  systematic  block  codes  in  the  Hamming  metric,  in 
particular  the  MCJ  codes  [1] .  The  MCJ  codes  of  immediate  interest  are  of  code  length  a 
prime  number  p.  Most  prime  number  alphabet  sizes  commensurate  with  record  communica¬ 
tion  requirements  and  with  reliable  computer  to  computer  communication  requirements 
appear  to  be  satisfied  by  primes  p  <  104 . 

The  procedures  discussed  here  are  quite  general,  but  the  description  is  given  in  terms 
of  a  specific  example  using  the  prime  number  p  -  37  in  the  finite  field  GF  (37).  This  is  the 
smallest  prime  number  that  can  contain  a  full  numeralphabet  including  the  26  letters,  10 
numerals,  and  a  spacing  symbol.  These  MCJ  prime  number  length  codes  have  superior 
error-correcting  properties  because  they  are  MDS  codes,  can  be  decoded  in  a  straightfor¬ 
ward  manner  using  well-defined  algebraic  techniques  including  the  theory  of  algebraic  and 
geometric  invariants  not  in  general  available  to  Bose-Chandhuri-Hocqueuhem  (BCH)  code 
structures,  and  can  be  chosen  to  closely  match  communication  circuit  and  user  requirements 
over  real  channels. 

As  BCH  codes  are  defined  by  generating  polynomials  with  distinct  roots,  and  as  the 
most  useful  forms  of  MCJ  codes  are  generated  by  a  repeated  root,  the  usual  BCH  synthesis 
techniques  prove  to  be  inadequate  for  these  codes.  However,  many  decoding  techniques 
usually  associated  with  BCH  codes  can  be  used  to  decode  MCJ  codes  [2-4] . 

Note  especially  that  the  MCJ  codes  imply  the  removal  of  the  usual  restriction  on 
cyclic  codes  that  the  code  length  n  and  the  characteristic  p  be  relatively  prime,  usually 
indicated  by  (n,  q)  =  1,  where  q  =  pm .  This  gives  the  MCJ  codes  a  slight  coding  advantage 
over  the  other  optimal  or  MDS  codes,  the  Reed-Solomon  (RS)  codes. 

To  demonstrate  the  general  principles  involved,  the  code  to  be  discussed  has  generating 
polynomial  g(x)  =  (x  -  c)n  "  k  and  will  be  described  in  the  cyclic  case  where  the  constant 
c  *  1,  the  code  length  n  =  37,  and  the  number  of  information  symbols  k  =  31. 


Manuscript  submitted  January  11, 1980. 
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The  starting  point  here  is  with  the  Vandermonde  matrix,  which  in  the  multiple  root 
case  is  a  generalization  of  the  Pascal  triangle.  This  matrix  V(c)  for  c  =  1  gives  the  powers 
(x  -  iy  ,  0  <  i  <  36.  Figure  1  emphasizes  its  relationship  to  the  system  function  S0,  as  it 
evolves. 

Figure  2  shows  the  location  of  the  basis  vectors  for  the  row  space  that  defines  all  the 
code  vectors  of  the  (37,  31)  code.  The  code  generator  matrix  G1(x)  is  developed  from 
multiples  (x  -  c)'  g(x)  instead  of  x‘g(x)  as  in  the  usual  synthesis  of  cyclic  codes.  To  be  noted 
specifically  in  Fig.  2  are  powers  (x  -  1)',  31  <  i  <  36,  which  when  read  as  vertical  vectors 
from  left  to  right  also  are  the  truncated  powers  of  (x  +  iy,  0  <;  <  36,  and  which  are  in  the 
form  of  a  matrix  HUx)  whose  rows  are  generated  by  (x  -  l)‘h(x),  and  h(x)  ■  (x  -  1  )nlg(x)  = 
(x-  1)31. 

The  matrix  in  Fig.  3  contains  the  reduced  echelon  [7feP]  as  derived  by  elementary 
row  operations  on  Gx  of  Fig.  2.  The  transmit  system  matrix  So  *  Sj  is  shown  in  Fig.  4, 
which  defines  the  systematic  MCJ  code.  The  decoding  receive  system  matrix  S~^  shown  in 
Fig.  5  satisfies  SqSq1  =  /,  where  I  is  the  37  X  37  identity  matrix.  The  matrix  generates  syn¬ 
dromes  exactly  as  in  Fig.  2  of  reference  [1] ,  but  physically  in  an  entirely  different  manner, 
even  though  both  are  implementations  of  (x  +  iy.  It  is  interesting  to  note  that  the  physical 
realization  of  the  preceding  encoders  and  decoders,  when  implemented  by  linear  (or  non¬ 
linear)  sequence  generators,  is  more  naturally  described  by  the  canonical  form  given  by  Bose 
and  Chaudhuri  in  their  original  paper  [3] ,  than  by  Peterson  and  Weldon  [4] .  However,  the 
Massey  virtual  encoder  (1]  is  best  described  by  Figs.  4  and  5  of  this  report. 

Several  different  equivalent  definitions  for  MCJ  codes  are  possible  as  approached  by 
different  points  of  view,  but  no  attempt  is  made  here  to  exhaust  these  possibilities.  Finally, 
although  the  MCJ  codes  will  effectively  correct  random  errors,  they  also  have  inherent  burst 
error  correction  capabilities. 

An  important  consideration  in  the  synthesis  of  block  codes  is  the  inclusion  of  fail-safe 
error  correction  procedures  in  the  decoding  algorithm  so  that  catastrophic  error  correction 
or  error  propagation  is  minimized.  The  following  theorem  is  included  as  one  way  of  accom¬ 
plishing  this;  however,  the  most  productive  test  is  outlined  in  paragraph  4.0  of  the  decoding 
algorithm. 

Theorem:  If  the  homogeneous  matrix  of  coefficients  is  of  rank  r,  the  nonhomogeneous 
matrix  or  system  of  coefficients  is  consistent  provided  that  the  rank  of  this  augmented 
matrix  is  also  r. 

Consistency  tests  in  the  decoding  procedure  are  pursued  under  five  hypotheses  as 
follows: 

H4 :  Four  or  more  errors 

tfg :  Three  errors 

H2 :  Two  errors 

Hx:  One  error 

H0:  Zero  errors. 
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Fig.  2  — Baai*  factors  for  (37, 31)  code 


Pig.  8  —  Reduced  echelon  matrix 
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The  specific  purpose  here  [2]  is  to  distinguish  between  cases  where  <  t  errors  occur, 
and  cases  where  >  t  errors  occur  and  correction  should  not  be  attempted. 

CORRECTION  PROCEDURE 
Preliminary  definitions: 

Let 

d  “n-*  +  l»2f  +  l 
m0  “0 
m  -  1 

n  *  pm  *  p  m  code  length 

k  ■  number  of  information  digits 
n-  k  *  number  of  parity  digits 
t  ■  maximum  number  of  correctable  errors 
Transmitted  signal 

fix)  -  f0+fl*  +  -+fn-k-l*n~k'1+~+fn-l*n~l- 
Noise  signal 

«(*)  *  eo  +  elx  +  •"  +  en-k  -  1** ”  '  1  +  —  +en-  ix"  ~  *• 

Received  signal 

r(x)  =  f(x)  +  e(x); 

r(x)  *  rQ  +  ri*  +  •••  +  rn-*- 1  ♦  ...  *r  1xn~ l. 

The  syndrome 

s(x)  «  *o  +glx  +  -+tn-k-l*n~k~1’> 
g(x)  *(x-c)p~k. 

Specific  Example 

The  following  procedure  is  adequate  for  the  decoding  of  any  MCJ  code  described  here 
and  uses  a  specific  example  in  the  Galois  field  GF(p)  with  p  =  37,  and  the  following  param- 
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eters:  c  =  1,  n  =  37,  k  =  31,  n  -  k  =  7,  and  t  =  3,  which  will  correct  up  to  three  character 
errors.  The  block  diagram  to  be  followed  is  given  in  Fig.  6. 

The  received  signal  sequence  is  entered  into  the  box  r(x),  following  which  the 
syndromes  si  are  calculated  in  the  syndrome  former,  more  explicitly  detailed  in  Fig.  7. 
The  modified  syndromes  Sg  in  Fig.  6  are  calculated  [5]  as  shown  in  Fig.  8  as  a  matrix 
multiplication  derived  as  follows:* 

St  =  13  |  .  J  >!c'V  1  <  i  <p  -  k 


•An  error  in  Eq.  (32)  in  Ref.  1  is  corrected  here. 
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Fig.  7  —  Syndrome  former 


(*0*1*2*3*4*5* 


'  1  0  0  0  0  0 

0  11111 
0  0  2  6  u  30 

0  0  0  6  36  2 

0  0  0  0  24  16 

0  0  0  0  0  9 


*  <S0S1S2S3S4S5) 


Fig.  8  —  Matrix  multiplication  for  deriving 
modified  syndrome* 
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For  m0  *  0,  and  noting  that  0°  =  1,  the  various  matrices  are  related  by 

Nt_j-Vt_JPt_jvlJ, 

where  the  distinct  root  Vandermonde  matrix  Vt_j  is  defined  as  follows: 


DECODING  ALGORITHM 

The  procedures  in  the  following  decoding  algorithm  are  summarized  in  Fig.  9. 
Step  1: 

1.0  Define: 
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2.1  If  B0  *  B1  *  B2  ■  0,  go  to  step  3. 

2.2  If  B0  =  0  and  at  least  one  Bi  0,  i  =  1,  2;  inconsistent,  do  not  correct. 

2.3  If  Bq  #  0,  then  find 

S2  fli 


Step  3: 

3.0  Define: 


Cl-|-S1l;C0-|50|. 


3.1  If  Cq  -  Cj  =  0;  no  errors. 

3.2  If  CQ  =  0  and  C1  ¥=  0;  inconsistent,  do  not  correct. 

3.3  If  CQ  #  0,  then  find 


-o. 


■;and  -ox  +* 


0. 


Step  4: 

4.0  Substitute  b  in  the  appropriate  equation,  0  <  b  <p  -  1,  for  x,  as  shown  in  Fig.  10  (the 
well-known  Chien  search),  and  find  the  error  locations x0,  Xj,  x2, ...» xf_  _  j.  Go  to 
step  5.  If  the  number  of  error  locations  does  not  equal  the  degree  t  -  j  of  the  proper 
equation,  more  than  t  -  j  errors  may  have  occurred;  inconsistent,  do  not  correct. 

4.1  Alternate  Step  4:  Solve  appropriate  equations  by  algebraic  techniques,  indicated  in  Fig. 
11,  to  find  error  locations. 


Step  6: 


5.0  Define: 


NRLRVOftTMM 


r\§.  11  —  Algebraic  technique  for  find  inf  error  locations 
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The  error  values  m: 


5.1  Define: 


The  error  values  are: 


NRL  REPORT  8393 


5.2  Define: 


IV,  I  -  111. 


Hie  error  value  is: 


y o  “ 


is„i 

iv,  i 


Step  6: 

6.0  As  a  final  step,  the  fully  corrected  text,  together  with  the  corrected  parities,  is  once 
more  entered  into  the  syndrome  former  for  final  verification.  The  new  syndrome 
will  now  be  zero,  indicating  either  the  message  has  been  successfully  corrected,  or 
mucorrected  to  another  codeword.  All  other  errors  at  distances  greater  than  d  will  have 
been  detected  in  previous  inconsistency  checks.  The  complete  procedure  is  outlined  in 
Fig.  6. 


SIMPLIFIED  8YNTHESIS 

In  practical  code  synthesis,  the  steps  shown  in  Figs.  1,  2,  and  3  may  be  bypassed  by 
first  defining  Sq1  ,  which  is  always  in  the  simple  binomial  coefficient  form  of  Fig.  5,  and 
then  finding  its  inverse  system  matrix  S0,  as  in  Fig.  4,  by  any  conventional  means.  Further, 
the  system  matrix  S0  may  be  independently  found  by  noting  that  most  of  the  parity  entries 
are  sequential  states  of  a  simple  nonlinear  generator  whose  coefficients  are  given  by  the 
expansion  of  (x  -  l)6;  and  the  initial  segment  given  directly  by  the  expansion  of  (x  -  l)1, 
i-0,1 . 6. 


CONCLUSIONS 

Several  alternatives  have  been  considered  for  software  or  hardware  implementation  of 
MCJ  codes  for  communication  system  error  protection.  It  was  emphasized  that  MCJ  codes 
provide  improved  flexibility  over  existing  procedures  in  matching  an  error  control  code  to 
a  given  user  alphabet  size  while  maintaining  optimal  performance.  Techniques  will  be  pre¬ 
sented  in  the  near  future, similar  to  those  given  here,  for  synthesis  and  decoding  of  Reed- 
Solomon  type  codes  also  over  the  Galois  ground  field. 
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